Reaction-Diffusion-Branching Models of Stock Price Fluctuations 
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Several models of stock trading [P. Bak et al, Physica A 246, 430 (1997)] are analyzed in 
analogy with one-dimensional, two-species reaction-diffusion-branching processes. Using heuristic 
and scaling arguments, we show that the short-time market price variation is subdiffusive with a 
Hurst exponent H = 1/4. Biased diffusion towards the market price and blind-eyed copying lead 
to crossovers to the empirically observed random- walk behavior (H — 1/2) at long times. The 
calculated crossover forms and diffusion constants are shown to agree well with simulation data. 
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The movement of stock prices is among the oldest class 
of fluctuation phenomena that have been analyzed quan- 
titatively , yet very limited understanding on the ori- 
gin and size of market volatilities — so fundamental to 
much of the finance literature || — is available to date. 
Recent empirical studies H|]|5| of major financial indices 
around the world have revealed a number of universal fea- 
tures in the time series, including Levy-like distribution 
of short-time price moves. In search of unifying princi- 
ples that explain the observed behavior, Bak, Paczuski, 
and Shubik (BPS) Q recently proposed several models 
of stock trading among interacting agents. Their numeri- 
cal study of these models suggests a rich set of dynamical 
behavior, including stock price fluctuations which exhibit 
similar statistical patterns as those of real markets || . 

From a statistical mechanical point of view, the family 
of models proposed by BPS are particularly interesting 
due to their connection to reaction-diffusion processes fa- 
miliar in physical contexts , and thus one may hope to 
gain some insight into the collective behavior of traders 
by exploiting this analogy. A mapping between the two 
is easily constructed by taking the target price of agents 
as their coordinates on a one-dimensional price axis. The 
two types of agents, i.e., buyers and sellers of a stock, are 
identified as two species, A and B, respectively. At any 
given moment, the population of buyers are separated 
from that of sellers by the market price xm(£) where 
transactions, or for that matter "reactions" A + B — > 0, 
take place. Such reaction-diffusion problems have been 
studied extensively in the past A result of partic- 

ular interest is the power-law scaling of the reaction front 
fluctuations, 



([x M (t)-XM(t')} 2 ) 1/2 ~\t-t'\ 
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where the exponent H = 1/4 (with possible logarithmic 
correction). There are, however, two new elements in the 
BPS models which have not been examined before: bi- 
ased diffusion of A and B particles towards the reaction 
front, and price "copying" which translates to branching 
A — > 2A and B — > 2B. BPS showed numerically that, in 



the latter case, the long time behavior of Xm(*) changes 
to that of a random walk with H = 1/2. 

The purpose of this paper is to establish an analytic 
foundation for various observations made by BPS in their 
pioneering work and also to further quantify and extend 
their numerical results. We identify the driving force of 
the market price variation and determine the size of the 
market response from the distribution of agents near the 
market price. The analysis yields not only the scaling 
exponent H, but also the scaling amplitudes and various 
crossovers. Good agreement is reached between theoret- 
ical predictions and simulation data for a broad range of 
model parameters. 

The original BPS model is a trading game with equal 
number (N/2) of buyers and sellers, each attaches a price 
Xi to the stock they intend to buy or sell. A buyer owns 
no share and a seller owns exactly one share. The agents 
perform simultaneous, independent random walks on the 
price axis until they meet a member of the other group. 
Upon transaction, buyer and seller exchange their role 
and are then relocated on the price axis. In this paper 
we shall consider three variants of the model as detailed 
below: 

Model I (unbiased diffusion) — The price Xi of agent i 
moves up or down by one unit with equal probability in 
each time step. For convenience, we decouple the reac- 
tion event A + B — > from the relocation of the agents 
in a transaction (see note p2| ). A steady-state situation 
is maintained by injecting new agents from the two ends 
of a prescribed price interval at a given rate J. 

Model II (biased diffusion) — The rules in this case 
are similar to those of Model I except that the updat- 
ing of Xi is biased towards the market price. Specifically, 
for a buyer, Xi — > Xi + 1 with probability (1 + D)/2 
and Xi — > Xi — 1 with probability (1 — D)/2. The rule 
is reversed in the case of sellers. Obviously, Model I is 
regained by setting D = 0. 

Model III (biased diffusion with copying) — The up- 
dating of Xi is the same as in Model II but now, after 
a transaction, the buyer and seller are immediately re- 
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injected into the market by duplicating the price of a fel- 
low agent chosen at random. According to BPS ||, such 
a process imitates herding behavior in real markets. In 
the particle language, it can be represented by stochastic 
branching A — > 2A and B — > IB. 

We start our discussion by considering a modification 
of the above models which is minor from the point of 
view of a given diffusing particle in its whole lifetime 
(i.e., from its first release to the reaction), but it triv- 
ializes the problem completely. Instead of asking an A 
particle to find a B particle for a reaction, we assume 
that the reaction always takes place at a fixed position, 
say x = 0. In essence, we are making the assumption 
that the market price fluctuates at a much slower rate 
compared to the diffusive motion of individual agents, a 
commonly used approximation in the study of interface 
fluctuations |l3[ . The point x = now serves as a trap 
of the diffusing particles which do not interact with each 
other. In the continuum limit, the average density of 
buyers a(x,t) and sellers b(x,t) obey the following linear 
equations, 



d t a = jd 2 a - (3d x a + S A , 
d t b = 7 dlb + (3d x b + S B , 



(2a) 
(2b) 



with the boundary condition a(Q,t) — 6(0, t) — 0. Here 
7 is the diffusion coefficient of individual particles and (3 
describes drift towards the current market price. The 
updating rule of Xi given above specifies f3 = D and 
7 = (1 — D 2 )/2 (parallel updating) or 7 = 1/2 (ran- 
dom sequential updating) . The source terms S A and Sb 
correspond to injection of new particles into the system. 
We now determine the steady-state solutions ao(x) and 
bo(x) = a (—x) to Eqs. (g). 

Models I and II. — The particle current 



(3) 



is a constant in this case. For (3 — (Model I), the 
current is maintained by a linear profile [Fig. 1(a)], 



a (x) = -Jx/j. 



(4) 



For (3 > (Model II), ao(x) crosses over from the linear 
function dJ) close to the origin to a constant J / (3 at large 
distances [Fig. 1(b)], 



a (x) = (J//3)[l-exp(/3x/7)]. 



(5) 



Model III. — Branching introduces source terms S A = 
aa and Sb — cub, where a is a branching rate. The equa- 
tion for ao (x) is now a second-order ordinary differential 
equation which can be solved to yield, 



ao(x) = C[exp(fc_x) — exp(fc+a;)] 



(6) 



where k± = ((3 ± y/ '(3 2 — 4aj) / (27) and C is an overall 
amplitude. The steady-state solution exists only when 



a < a c = (3 2 /A 7 . The shape of the profile is indi- 
cated in Fig. 1(c), which is linear close to the ori- 
gin and decays exponentially at large distances. The 
current of incoming particles at the origin is given by 
J = -jd x a = -yC(k + - fe_). 
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(a) (b) (c) 

FIG. 1. Average density profiles of buyers (do) and sellers 
(bo) along the price axis for (a) Model I (unbiased diffusion); 
(b) Model II (biased diffusion); (c) Model III (biased diffusion 
with copying). Note that in all three cases ao(x) and 60 (a;) 
vanish linearly at the market price xm- 

The total number of, say A particles that arrive at 
the trap in a time interval t to t' = t + r is a fluctuating 
quantity which can be expressed as, 



n A 



(7) 



Here the sum runs over all particles i that entered the 
system from the beginning of the process to time t' , and 
rji is a random variable which takes the value one if parti- 
cle i is trapped during the interval r and zero otherwise. 
Since the rji's are independent from each other fl4|| , we 
easily find, 



(n A ) = 



(n|> - (n A f = 



Pi) 



(8a) 
(8b) 



where ^ is the probability that rji — 1. ¥ov pi <C 1, which 
holds when r is much smaller than the typical spread of 
the lifetime of the diffusing particles, we have the follow- 
ing approximate relation, 



(n A ) - (n A ) 2 ~ (n A ) = Jr. 



(9) 



where, as before, J is the flux of particles entering the 
trap. Results derived below are based on this approxi- 
mation but other situations may also be considered. 

We now construct a heuristic argument to show how 
the fluctuations in n A and ns lead to a shift of the re- 
action front or the market price. To be definite, we take 
An = n A — us > so that an upward move of the mar- 
ket price from xm at time t to x' M at time t' = t + r 
is expected (see Fig. 2). Loosely speaking, the interval 
[xm,x' m ] defines a reaction zone within which most of 
the n A particles entered through xm reacted with the 
tlb particles entered through x' M . The excess number of 
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A particles (buyers) An have either reacted with the B 
particles initially in the zone at time t, or remained in the 
zone at the end of the period. Based on this observation, 
we may identify An with the sum of shaded areas under 
b(x,t) and a(x,t') in Fig. 2, respectively, 

r x M r x M 
An ~ / a{x,t')dx + / b(x,t)dx 

Jx M Jxm 
i-Ax 

~ 2 / a {x - Ax)dx. (10) 
Jo 



Here Ax = x' M — xm is the price move over the time 
period r. A similar relation holds for An < 0. 

n A n B 




%m %m price 
FIG. 2. Schematic illustration of how an excess number of 
particles An = tia ~ ns reaching the reaction zone [xm,x' m ] 
is accommodated by a shift of the density profiles at time 
t to those at a later time t'. Graphically, An is identified 
with the sum of the shaded areas under a(x,t') and b(x,t), 
respectively. 

Equation JTo|), which holds in an average sense, is our 
fundamental relation that links the price move Ax to the 
"demand-supply imbalance" An through the density pro- 
file ao(x). From the statistics of An we can then work 
out the statistics of Ax. Below we give our results on the 
market price fluctuations using this approach. 

Model I. — In this case, the intrinsic profile oq(x) is lin- 
ear. From Eqs. (0) and (|lC|), we obtain An ~ J(Ax) 2 /7. 
From the variance (An 2 ) = 2Jt, we obtain, 

([x M (t)-x M (t')]y /2 ^(^) 1/4 \t~tf/\ (ii) 

This form agrees with previous results of Refs. [| -0. 
(At very short times, the discreteness of An manifests 
itself which leads to deviations from Eq. (|TT|) . See also 
Ref. §.) 

Model II. — Since ao(x) goes to a constant for \x\ > 
x c = 7//?, Eq. (|l0|) yields a linear dependence An ~ 
2JAx/(3 for Ax > x c . Hence the move of x^if] at 
long times is a random walk with a diffusion constant 
r ~ (Ax 2 )/2r = /3 2 /(4J). More detailed calculation 
yields a crossover scaling, 

([x M (t) - x M (i')] 2 > 1/2 * x c <J>(^% (12) 



where t c — 7 2 J//3 4 . The limiting forms of the scaling 
function are given by <&(s) ~ (As/n) 1 / 4 for s C 1 and 
$(s) ~ (s/2) 1 / 2 for s > 1. 

Model III. — In this case the profile (||) extends only 
over a finite range of x, so the finite lifetime of a par- 
ticle becomes an important factor in our consideration. 
The short time behavior of the price fluctuation is sim- 
ilar to that of model I and II due to the linear behav- 
ior of ao(x) close to the origin, which is common in all 
three cases. Thus Eq. ( |ll"| ) can still be applied in this 
regime. Crossover to a different behavior is expected 
when t = t' — t becomes comparable to the lifetime of a 
particle To = 2j/f3 2 . In fact, To is also the relaxation time 
of the density profiles as can be seen by bringing Eq. (^) 
into a dimensionless form. On time intervals larger than 
To, memory about the initial profile is essentially lost and 
the next move of the market price is equally likely to be 
up or down, hence a random walk behavior with a step 
size set by the size of the fluctuation x = 7 1 / 2 J -1 / 4 ^ 4 
at t = To [see Eq. (Ell)]. The usual scaling argument 
then yields, 

a*MW-*M(to]y /2 - (13) 

V T / 

where *(s) ~ s 1/4 for s < 1 and *(s) — s 1/2 for s > 1. 
The diffusion constant of the market price at long times 
is given by V ~ x^/t = /?7 1/2 J~ 1/2 . 

We have performed numerical simulations of the BPS 
models to check the validity of the theoretical analysis 
presented above. Since our results on Model I are sim- 
ilar to those of previous studies (apart from a possible 
logarithmic correction), we shall focus on Model II and 
III. 

Model II was simulated at J = f3 where the asymptotic 
density is one as in Ref. @. The system size is chosen to 
be N = 2000 or larger to ensure that the reaction front 
does not fluctuate out of the boundaries during the time 
period simulated. Otherwise, N is found not to have any 
significant effect on our results The system is first 
equilibrated for a period to = 10to where To = N/(2(3) is 
the typical lifetime of a particle. The market price time- 
series xjvf(i) is then recorded over 500 successive time 
segments, each of length 8192 time steps. We then cal- 
culate {[xM(t) — xm(£')] 2 ) averaged first over each time 
segment and then over different segments. In Fig. 3(a) 
we plot the simulation results using scaled variables for 
D = 0.01 to 0.5. There is indeed a good data collapse 
over six decades. In fact, for t > t c , not only the scaling 
exponent, but also the scaling amplitude are borne out 
by the data. 

The simulation of Model III was carried out in a similar 
way as that of Model II, except the number of particles 
N is now fixed. To compare the simulation data with Eq. 
(|l3|), we use the relation J — aN/2 (particle conserva- 
tion) from the solution (J^). Taking a — a c = /3 2 /47 [ p"5| , 
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we obtain xq = 2 r )f3~ 1 N~ 1 / i . In Fig. 3(b) we plot the 
simulation results for market price fluctuations using the 
scaling suggested by Eq. (|f|). For four different values 
of D = P and two system sizes N — 400 and 1000, good 
data collapse is again achieved. 
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FIG. 3. Scaling plots of the market price fluctuations ver- 
sus time increments for a wide range of model parameters, 
(a) Model II. The value of D for each data set is given in the 
legend, (b) Model III. The values of D and N for each data 
set are given in the legend. 

In summary, we presented a heuristic method to link 
the market price fluctuation to the diffusive motion of in- 
dividual agents using the BPS models as examples. The 
analysis yields qualitative as well as quantitative predic- 
tions on the size of the market price fluctuations as a 
function of time, the number of traders in the market, 
and various other model parameters. For short times, a 
previously known H = 1/4 scaling law is rederived and its 
validity is correlated to the generic linear shape of pop- 
ulation density profiles near the market price. Crossover 
to the long-time random walk behavior with H = 1/2 
takes place when agents are driven to the market price 
via a diffusion bias. Expressions for the crossover time 
as a function of various model parameters are derived. 
These results are shown to compare favorably with the 
simulation data. 

The H = 1/4 scaling at short times is quite remark- 
able and is against the prevailing thinking in finance that, 
in a market with noise traders only, there should be no 
restoring force to price moves and hence no correlation in 



the market price time series. Although the exponent is 
not new, the analysis presented here makes it plain that 
resistance to price change is inherent in the existing price 
distribution of agents. It remains to be elucidated how 
such tendencies are modified when external information 
(e.g., financial news) are fed into the market. 
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Foundation for financial support and the Physics Depart- 
ment, Hong Kong Baptist University for their hospitality. 
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